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I. INTRODUCTION 



In recent years the great expansion in the development 
of semiconducting materials with acceptable thermoelectric 
characteristics has made them attractive for, in special 
applications, direct energy conversion and heat pumping. 

Thermoelectric power sources and heat pumping devices 
can offer significant advantages with respect to long life, 
years of unattended operation, freedom from environmental 
restraints, and unusually favorable specific power ratios. 

It is true that even now thermoelectric cooling and 
power sources are not serious competitors of powerful heat 
pumps and conventional power apparatus, but they are 
irreplaceable in those cases where the overriding require- 
ments are small size, simplicity of construction, quiet 
operation, weight, etc. 

Various microrefrigerators, used in radio engineering, 
medicine and industrial design, have received proper 
recognition. Fueled portable thermoelectric generators, 
especially for military applications, are already used. 

In the present work an attempt is made to design, at 
least theoretically, a thermoelectric air cooling device 
servicing the passenger compartment of a small automotive 
vehicle, in conjunction with a thermoelectric generator 
using the waste heat of the exhaust gasses of the engine 
to provide a portion of the required electrical power. 
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The feasibility of the replacement of the conventional 
alternator by a thermoelectric generator is also examined. 

The introduction of the thermoelectric generator will 
improve the coefficient of performance of the thermoelectric 
cooler, and will also increase the initial installation cost. 

The design was developed from the macroscopic point of 
view, considering the present state of the art of thermo- 
electric materials. 

The total work is divided into five sections: 

1) The heat flow analysis for the semiconductor material 
under a temperature gradient and internally generated heat 
due to the flow of electric current. The results of this 
analysis are almost common to the thermoelectric cooler and 
generator. In this section a differential equation, which 
describes the heat balance in the thermoelement has been 
derived. In this equation, the cross section of the thermo- 
element is treated as a function of the length and the 
properties of the material, as thermal conductivity, Thomson 
coefficient, etc., are treated as functions of temperature. 

2) The analysis and design of the thermoelectric cooler. 
A detailed analysis of the T.E . cooler is made for both modes 
of operation, maximum heat pumping rate and maximum coeffi- 
cient of performance. The analysis is closed with comparison 
between the two modes and the selection of max COP as the 
most suitable for the present application. 

The analysis is followed by an indicated design of a 
cooler with a capacity of 7000 BTU/hour which, from a quick 
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estimation, was found to be satisfactory for the selected 
environmental conditions. 

3) The analysis and design of the thermoelectric 
generator. The analysis is conducted from the standpoint 
of maximum efficiency mode. It is followed by selection of 
the proper thermoelectric materials for the assumed conditions . 
Computer techniques are employed for estimation of the 
temperature relations of the thermoelectric material proper- 
ties. The analysis is closed with the estimation of the 
available heat power from the exhaust gasses. The attempted 
design concerns a segmented type thermoelectric generator 
which is carried out exclusively by a computer program. 

4) Conclusions and recommendations for further study. 

5) Appendices 

Generally, the proposed configuration would be as it is 
shown in Figure 1. 
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Control Device 




FIGURE I 

DIAGRAMMATIC ARRANGEMENT OF THE PROPOSED SYSTEM 
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II. STEADY STATE HEAT FLOW IN THE SEMICONDUCTOR ELEMENT 



The derivation of the heat flow balance equation through 
the thermoelement is presented when the cross section of the 
element is a function of the length and the characteristics 
of the material, electrical resistivity and thermal conduc- 
tivity are functions of temperature. 

Under the proposed geometry of the semiconductor, the 
one dimension heat flow is evident. 




FIGURE 2 

PROPOSED GEOMETRY OF THERMOELECTRIC SEGMENT 
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The semiconductor is considered to be adiabatically 
insulated on all sides except side ABCD and EFGH . 

The temperature of the solid is assumed to be T__ at the 

H 

front side (ABCD) and T_ at the rear side (EFGH) , when T„ 

v~ n 

is greater than T c » The characteristics of the material 
are considered to be variable and more specifically to be 
functions of temperature which can be expressed as: 

1) Conductivity ^\(T) = a Q + a^ + a 2 T 2 

2) Resistivity p (T) = b Q + b 1 T + b 2 T 2 

o 

3) Seebeck coefficient o( (T) = c Q + c^^T + c 2 T 

da ( T} 2 

4) Thomson coefficient 6“(T) = T ■ ■ = c^T + 2c 2 T 

Considering the part of the semiconductor which is 
confined between the cylindrical surfaces of radii r and 
r+dr, let q(r) and q(r+dr) be the rates of heat flowing in 
and out of the differential element. 

If a current I is flowing through the semiconductor 
Joule heat will be generated due to the electrical resistance 
of the semiconductor and heat will be released or absorbed — 
depending on the direction of the flowing current — due to 
the Thomson effect. The differential of the Joule heat 
generation rate is given as : 




( 1 - 1 ) 
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where A(r) is the surface of the differential solid element 
at radius r. 

Also the differential form of the Thomson effect is 
given as : 



^ - I&(T)d7 (1-2) 

Consider now, for balance the rate of heat flow out of 
the differential solid element must be equal to the rate of 
heat flow in, plus the rate of heat generation in the element. 
These can be expressed as: 

+ ^ d-3) 

But q(r+dr) can be expanded in a series of the form: 

+ + . .. (i_4) 

If the first two terms will be considered, then equation 
(1-3) can be written as: 

^( r ) clr = q ± <) (1-5) 

>T 

Substituting the equivalent relations for qj and q T equation 
(1-5) becomes 

= + jg-(T) il— d-6) 

dr A(<”) dT " 
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The conduction heat — rate of flow — q(r) is given by: 



<j,(r) = - y(J)A(r) 



JT 

dr 



(i 



By differentiation in respect to r equation (1-7) becomes: 



JlACTfljT A/^clT v ,JLA0r)]dT 

dr ~ dT dr m ; dr ~* a;- 3F'dr'“ ^ ( 



dr 



dr* 



By substituting of equation (1-8), equation (1-6) becomes: 






For the geometric configuration of Figure 2 the cross 
section area A(r) is a linear function of radius r: 

ACO - (i 



ACO = nu- 



ll 



where m = 



Using the above expression for the cross section, 
equation (1-9) becomes: 



(1 



-7) 

- 8 ) 

-9) 

- 10 ) 

-10a) 

-ID 
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Equation (1-11) is the differential equation which 
describes the rate of heat flow through the semiconductor 
material, with the geometry of Figure 2 and the conditions 
stated before, when a temperature gradient is impressed and 
an electrical current is allowed to flow. 

This equation is non-linear and at least a rigorous 
analytic solution is not obvious. 

If in the range of temperature, when the thermoelectric 
couple is operated, the characteristics of the semiconductor 
material are constant or reasonably can be assumed to be 
constant, then equation (1-11) is greatly simplified and 
takes the following form: 



Jjl + _l + _ o 

dt** r dr 



( 1 - 12 ) 



Since the temperature gradient under which the thermo- 
electric element of a cooler is operating is small, equation 
(1-12) reasonably can be used as the basic equation for the 
design. 

This is not the case for the thermoelectric generator 
where the range of temperature is quite extended, and the 
material characteristics present sufficient deviations from 
the corresponding values at low temperatures. 

For this, equation (1-11) will be used as the basis for 
the thermoelectric generator design. 
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An analytic solution of equation (1-11) is given in 
APPENDIX A under certain assumptions. APPENDIX B provides 
the analytic solution of equation (1-12) . 
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III. THERMOELECTRIC COOLER 



A. ANALYSIS OF THE ELEMENTARY COOLER 

The operation of the thermoelectric cooler is based on 
the Peltier effect, which was discovered in the eighteen- 
thirties. This phenomenon was explained correctly by Lenz 
in 1838. When an electric current I flows through a contact 
between two conductors, a certain amount of heat is either 
evolved or absorbed at the contact, depending on the 
directions of the current. This heat is known as the Peltier 
heat Q n • 

Becquerel and others established that Q n is directly 
proportional to the current I : 

Q. n =TKY)I (3-D 

where TT (Y) is the so-called Peltier coefficient. 

It has also been shown that if the current flowing 
through the junction has the same direction as heat flow 
which is due to the heating of this junction then the 
Peltier heat is absorbed; while in the opposite case, heat 
is evolved. 

The Peltier effect has been observed in a great variety 
of materials, and the Peltier coefficient has been found to 
be dependent on temperature in a complex fashion. 
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The thermocouple is considered under steady state condi- 
tions with its hot junction at a temperature T H and the cold 
junction at 

Electric current is considered flowing through the couple 
in the direction opposite to its Seebeck voltage. 

The arrangement is shown in Figure 3. 

The law of conservation of energy gives: 









(3-2) 



where q Q is the rate of heat absorption at the cold junction, 
is the rate of heat which is evolved at the hot junction, 
and W is the electric power. 

An important efficiency index of the thermocouple is the 
coefficient of performance which is defined by: 



COP - 




(3-3) 



Another important index is the heat pumping rate — or 
the refrigerating effect — q Q which can be expressed by: 



\=TC)I + Y T) A p (r) 57 + X(t) A*Cr) 57 - 



(3-4) 



a 



The thermoelectric coolers are usually operating under 
a small temperature difference. This is particularly true in 
the air conditioning applications. Under these small 
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FIGURE 3 

ELEMENTARY THERMOELECTRIC COUPLE 
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temperature ranges, the variations of the properties of the 
semiconducting materials are limited, and the use of average 
values introduces slight errors of the order of one percent. 

The use of average values simplifies the analysis of the 
system and they are used throughout this section. 

From APPENDIX B the rate of heat flowing to the cold 
junction, because of the existing temperature gradient and 
the generated Joule heat, for the P-type and N-type segments 
is: 






-\M) 

A /v A' J f- 



(±Ek AT-ifp 
r 5 ff 

_ a (4l)/v A T _ i. \ ^p 

‘& 1 ( 4 ^)„ 



(3-5) 



(3-6) 



Introducing the following variables: 



X = 




3 = 




(3-7) 



(3-8) 



n(T c >- ( U f (T c ) l 4 - [cUTc)l ) T c 
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equation (3-4) becomes: 



v- - ir[e p x + a] o-io 

The work done by the electric source, per unit time, 
consists of the Joule heat evolved in the thermocouple arms, 
and the work done against the thermo-e.m.f To(AT 

Ntf = IMT +I 2 K (3-11) 



where R is the ohmic resistance of the couple. The coeffi- 
cient of performance can be written as: 



iruvx+en] 
w&t .ne f x*fo] 

co?= 5 ^ 

I*A74l 2 K 



(3-12) 



(3-13) 



B. CONDITIONS FOR MAXIMUM HEAT PUMPING RATE 

The heat pumping rate of the elementary thermoelectric 
cooler, as it appears in equation (3-10) , is a function of 
X, Y, and I. Relative maximum can be found by the Lagrange 
multiplier technique using as constraint function: 

^=p P **e,rJ (3_14) 
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Thus the following conditions must be satisfied: 



4 U = o 

s>x ' p 



1% 






llo 

93 



= 0 



(3-15) 

(3-16) 

(3-17) 



where h is the Lagrange multiplier. 

Eliminating h from equations (3-15) and (3-16) the first 

condition for q to be maximum is found to be: 
o 




By substitution of equations (3-7) and (3-8) , equation (3-18) 
becomes : 



(^t)* - [_ft 



;,v X k 
e.%. J 



(3-19) 



The second condition for the maximum heat pumping rate is 
found from equation (3-17) : 



lo^L — 






(3-20) 



The voltage to be applied to the thermocouple is given by: 



V=* AT+io^ 



(3-21) 
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Substituting the expression for I from equation (3-2) in 
the above equation it is found: 



V= o( T h 



(3-22) 



The required electrical power for the couple is : 

p - VI - (X l H Io ft . (3-23) 



The maximum heat pumping rate q is found from equation 

max 

(3-10) after the substitution of equations (3-18) and (3-20) . 






non cxx 



o( Z Tc (P~J) (D~*J) 
2d (D-l)(D + 0 



(3-24) 



where 






(3-25) 



1 = 









Thh-Tc 






Th 

Tc 



(3-26) 



(3-27) 
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The coefficient of performance under the maximum heat pumping 
rate conditions is found from equation (3-13) after trans- 
forming it by means of equations (3-18) and (3-20) . 



CO-3HP+5) 



(3-28) 



23 (o-l)(D + l) 



The rate of heat rejection is determined with the help of 
equations (3-3), (3-28) and (3-24). 



C. CONDITIONS FOR MAXIMUM COEFFICIENT OF PERFORMANCE 

As it was in the case of maximum heat pumping rate, the 
coefficient of performance equation (3-13) , is also a function 
of x, y, and I. 

Applying the same technique, the following conditions 
must be satisfied: 




(3-29) 




(3-32) 



(3-30) 



(3-31) 
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Where the notation K instead of COP, for the coefficient of 
performance is introduced for convenience, and h again 
represents the Lagrange multiplier* Elimination of h between 
equations (3-30) and (3-31) gives exactly the condition of 
equations (3-18) and (3-19) • 

By algebraic manipulation the following expressions are 
obtained: 



t ck AT 


(3-33) 


(D-l)K 




w R At D 


(3-34) 


(D-i) 




r r\ -» 




Tc t ^ ~ Tc ] 


(3-35) 



It is obvious that equation (3-25) goes to zero when: 



D = — 

Tc 



(3-36) 



and it goes to infinity when: 



AT = 0 



(3-37) 



Combining equations (3-10), (3-18) and (3-33), the heat 
pumping rate becomes: 



V 






(3-38) 



R(D-a)* ( D *l) 
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D. COMPARISON OF THE EXTERNAL CONDITIONS 



The two cases of maximum heat pumping rate and maximum 
coefficient of performance differ very markedly as regards 
both the working conditions — voltage required — and the 
efficiency indices . 

Comparing the values of <3 max and q Q in Figure 4, it is 

seen that <3 max is greater than q Q and consequently, in order 

to have the same heat pumping rate for thermocouples working 

under conditions of q and COP it is necessary to use a 
max max 

larger number of thermoelectric units in the second case than 
in the first. This aspect is not decisive because from the 
efficiency point of view it is first of all necessary to 
ensure that the coefficient of performance is sufficiently 
high, since it is this coefficient which governs the elec- 
trical power consumed in a cooling device. On the other hand, 
it follows from Figure 5 that the coefficient of performance 

at the COP is greater than COP and that this difference 
max ^ 

can not be compensated at all. The preferable conditions 
are those when the coefficient of performance is a maximum 
since the materials known at present do not give such a high 
coefficient of performance that it can be afforded to have 
it lower, even for the purpose of economy of material. 

Designs based on the maximum heat pumping rate could 
be used only for the rapid cooling of an enclosure or for 
cooling it to a lower temperature at the cost of additional 
electrical power consumption. 
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FICURE 4 



COMPARISON OF q max - CURVE 2 - AND q Q - CURVE I 
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Coefficient of Performance 




FIGURE 5 

COMPARISON OF COPj^ - CURVE I - AND COP - CURVE 2 
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E. ESTIMATION OF TEMPERATURES OF OPERATION 

Complete air conditioning involves the simultaneous 
control of the following factors: 

1) Temperature 

2) Humidity 

3) Air motion 

4) Air distribution 

5) Dust 

6) Bacteria 

7) Odors 

8) Toxic gasses 

Of these the first three, temperature, humidity, and 
air motion, are paramount, and no air conditioning system 
deserves the name unless it satisfactorily controls those 
three factors. 

In the present work the thermoelectric air conditioning 
is examined only from the standpoint of the air cooling 
process . 

The device should be also considered as equally appli- 
cable for air heating process by just reversing the direction 
of the current through the thermoelements . 

From the standpoint of human comfort, the object of the 
air cooling device is not to cool the body but to control 
its rate of heat loss within reasonable limits by methods 
that are known to be the most conductive to comfort. 

There are four factors which, by their combined effect, 
determine the cooling power of the atmosphere upon the human 
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body. They are temperature, humidity, air motion and 
radiation . 

The effect of humidity and air movement are matters of 
common experience. When the indoor air is very dry, a 
higher temperature is necessary for comfort because of the 
greater cooling effect of the dry air. The evaporation of 
the perspiration with which the human body is always par- 
tially covered has a marked influence upon the rate of heat 
loss, and the rate of evaporation is dependent upon the 
humidity of the surrounding atmosphere. 

The cooling action of air motion upon the skin increases 
heat removal by the increased evaporative effect and by an 
increase in the convection effect as the layer of air in 
contact with the skin is repeatedly renewed. 

Loss of heat by radiation from the exposed body and 
clothing depends on the temperature of the surrounding walls 
and is independent of air temperature. 

Air movement reduces to some extent the radiation loss 
by lowering the surface temperature of the body. The rate 
of heat loss varies according to Steffen's fourth power law 
and is dependent upon the difference in temperature between 
the surface of the body and the mean temperature of the 
environment. This latter temperature is referred to as the 
mean radiant temperature. 

Experiments made on a large number of human subjects 
have demonstrated that there is a consistent relation between 
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the temperature, humidity and air motion in their effect 
upon comfort. When the results of such observations are 
plotted on a psychometric chart, it is found that the points 
representing conditions of equal warmth fall in approximately 
straight lines. These straight lines are called "effective 
temperature lines". For moving air the effective temperature 
for a given temperature and humidity is lower than that corres- 
ponding to still air. The effective temperature is the 
only true index of human comfort that is now available. 

Neither the dry-bulb nor the wet-bulb thermometer alone 
trully indicates the conditions of the air from a comfort 
standpoint. 

The optimum conditions for human comfort at summer time 
are considered to be 71 F° effective temperature and air 
movement of 15 to 35 feet per minute for an outdoor temperature 
of 80 F° . 

The following table gives the corresponding effective 
temperatures for several outdoor dry-bulb temperatures. 

OUTSIDE DRY-BULB EFFECTIVE 

TEMPERATURE (F°) TEMPERATURE (F°) 



105 

100 

95 

90 

85 

80 



75.5 



75 

74 

73 

72 

71 



TABLE 1 
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These effective temperatures are the statistical results 
of extensive experimental work, where several factors were 
considered, between them the stress which a human body suffers 
when it is moving from low to high temperature environments 
and vice versa. 

The outdoor conditions to be assumed in the calculations 
of the cooling and dehumidifying load should not be the 
highest recorded temperature in a location which usually 
exist for only a few hours during the year. They should be 
chosen as conditions which may be exceeded for a few days 
during the summer. 

In every large city a certain design dry-bulb and wet- 
bulb temperature which are readily ascertained, have become 
established and generally adopted. The following table gives 
these temperatures for a few cities of the U.S.A. 



LOCATION 


DESIGN 


TEMPERATURES 




DRY -BULB 


WET-BULB 


New York City 


95 


75 


Chicago 


95 


75 


Philadelphia 


95 


78 


Detroit 


95 


75 


Cleveland 


95 


75 


St. Louis 


95 


78 


Washington 


95 


78 


New Orleans 


95 


79 


San Francisco 


90 


65 


Dallas 


100 


78 


Jacksonville 


95 


78 



TABLE 2 



33 



Since an automobile is not designed to operate in a 
certain place, the worst case outdoor conditions are selected 
as those corresponding to the city of Dallas. 

So, outdoor dry-bulb temperature is 100 F° and indoor 
maintained-ef f ective-temperature the corresponding 75 F°. 
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IV. COOLER DESIGN 



A. GENERAL CONSIDERATIONS 

The design is based on the worst case which corresponds 
to an outdoor temperature of 100 F° and an effective tem- 
perature, for human comfort, of 75 F°. This means that the 
temperature of the incoming air in the hot side of the 
thermo elements would be 100 F° and the maintained tempera- 
ture in the passenger's compartment of the vehicle would 
be 75 F°. 

The proposed configuration of the thermoelectric cooler 
is shown in Figures 6a, 6b. It consists of radially 
arranged thermocouples, forming rings. A number of such 
rings are assembled to form the annular configuration. Fins 
have been introduced in both sides of the thermocouples to 
increase the rate of heat transfer. The outer annular finned 
space forms the heat absorbing exchanger and the inner cylin- 
drical finned space forms the heat rejecting exchanger. The 
cooling medium is air which flows by means of fans, one for 
every heat exchanger. When the vehicle is moving with a 
speed higher than a certain level, there is no need for 
operation of the heat rejecting exchanger fan. 

It is also assumed that the vehicle has been parked for 
a long period under the sun and the temperature of the indoor 
air has reached 100 F°. 

The temperatures of the hot and cold side of the thermo- 
elements, T^, T c , can be estimated if the desired logarithmic 
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FIGUEE 6a 

CROSS SECTION OF T.E. COOLER 
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FIGURE 6b 

LONGITUDINAL CROSS SECTION OF T.E. COOLER 
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mean and output temperatures of the cooling air, in both 
sides of the thermocouples are selected. 

According to the definition of mean logarithmic 
temperature : 




” ft 



'-iru -tl 




(4-1) 



the hot and cold shoe temperatures of the thermocouples are 
estimated as: 



Tc.rt- 



AT 0U ±_ T 



iotl 



A-i 



(4-2) 



where: 



A 

AT 

K 



AT 

= a k 



= T. - T 

in out 



= desired value of logarithmic mean temperature 



0 . . 0 , 



temperature differences as shown in 
Figure 7. 



The cooling capacity of the thermoelectric cooler is 
selected to be 7000 BTU per hour. This figure is considered 
as a moderate value, sufficient to serve a small vehicle 
carrying four passengers under the climatological conditions 
described previously. 
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ASSUMED TEMPERATURE DISTRIBUTION FOR 
DETERMINATION OF LOG. MEAN TEMPERATURES 
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As the most suitable thermoelectric material for cooling 
application selected a bismuth — tellurium — antimony com- 
pound for the P-type semiconductor, and a bismuth — tellutium — 
selenium compound for the N-type. These materials have been 
developed since 1968 from Westinghouse . The average proper- 
ties of the materials are displayed in Table 3. 



ITEM 


UNITS 


P-type 


N-type 


Seebeck Coeff. 


V/deg K 


205 x 10" 6 


- 220 x 10~ 6 


Resistivity 


ohm- cm 


1.2 x 10 -3 


1.08 x 10 -6 


Th. Conductiv. 


watts/deg 


c-cm 0.0121 


0.0149 


Fig. of merit 


1/deg K 


2.82 x 10 -3 


3.0 x 10 -3 



TABLE 3 

In order to obtain good mechanical strength the materials 
are processed by pressure and sintering. 

B. COOLER DESIGN BASED ON MAXIMUM COP 

This design is based on the results obtained from the 
analysis in Section III-C. 

1 . Thermoelectric Couple 

The coefficient of performance, as was defined in 
Section III-C, is independent of the geometric configuration 
of the thermoelectric couple. It is only a function of hot 
and cold shoe temperature and of the material properties. 

The heat pumping rate and the optimum current depend on the 
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total electrical resistance of the thermocouple which is 
a strong function of geometry. From this standpoint, an 
optimization of the geometry had been attempted using 
iterative computer techniques and the selected dimensions 
in Table 4 . 



ITEM 


UNITS 


P-type segment 


N-type segment 


R 1 


cm (in.) 


7.0 (2.756) 


7.0 


(2.756) 


r 2 


cm (in.) 


8.5 (3.346) 


8.5 


(3.346) 


4> 


degrees 


10.0 


o 

rH 

i— 1 




T 


cm (in.) 


1.0 (0.3937) 


1.0 


(0.3937) 






TABLE 4 







The calculated cooling and electrical 
of the thermocouple are shown in Table 5. 


characteristics 


ITEM 


UNIT 


VALUE 

-3 


Electrical Resistance 


ohm 


2.31 x 10 


Hot shoe temp. T^ 


degr. K 


322.73 


Cold shoe temp. T c 


degr. K 


291.5 


Electrical current 


amperes 


14.4 


Applied Voltage 


volts 


0.048 


COP 


dimensionless 


1.13534 


Required Power 


watts 


0.691 


Heat pumping rate 


watts 


0.74466 


Heat rejection rate 


watts 
TABLE 5 


1.43569 
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For the cooling load of 7000 BTU per hour a total 



number of 2754 thermocouples is required. If a nine degree 
wedge of insulation is provided for each thermocouple, 
then 12 couples form one ring and the total number of rings 
is determined as 230. Each ring, according to Table 4, 
has width 1.0 cm. Allowing insulation of 0.3 cm between 
the rings, the total length of the configuration is 299 cm. 



consisting of 57.5 rings, will satisfy the requirements. 
The externally supplied electrical power is calculated 
to be: 



and the heat rejection rate: 

2754 x 1.43569 = 3954 watts or 13504 BTU/hour 

2. Cold Side Heat Exchanger 

The introduction of extended surfaces, in the form 
of fins with rectangular cross section, increases the 
exchange surface of the heat exchanger. The fins have an 
efficiency given by: 



Four units, each having a length of 74.75 cm and 



2754 x 0.691 = 1903 watts 




(4-3) 
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where 



(4-4) 

Figure 8b shows the variation of n f as a function of 
for several values of h. The cooling medium is air with 
average properties displayed in Table 6. 



SYMBOL 


ITEM 


UNITS 


VALUE 




Thermal conductivity 


BTU/f t-hour-F° 


0.015 


P 


Dynamic viscosity 


lbm/f t-sec . 


1.7 x 10” 5 


P 


Density 


lbm/f t^ 


0.075 


C ? 


Thermal capacity 


BTU/lbm-F° 


0.24 



TABLE 6 



Figure 9 shows the heat exchanger arrangement. The long 
fins form four different types of air ducts, namely A, B, 
C, and D. 

Each air duct has an equivalent hydraulic diameter 
given by: 



D 



e 



4 S 
P 



(4-5) 



where S is the cross section and P the perimeter of the 
duct. 
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FIGURE 8a 



LONGITUDINAL FIN ARRANGEMENT 



Efficiency 
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FIGURE 8b 
FIN EFFICIENCY 




45 







FIGURE 9 

ELEMENTARY SECTION OF COLD SIDE HEAT EXCHANGER 
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The calculated equivalent diameters for each duct 



type are : 



(D ) , = 48.284 x 10“ 3 ft 

© A 

(D ) n = 43.89 x lO -3 ft 
e Jj 

(D e ) c =43.89 x 10" 3 ft 
(D e ) D = 62.26 x 10 _3 ft 



The selected temperatures for the air and thermo- 
couple cold shoe are displayed in Table 7. 



ITEM 

Incoming hot air 
Outgoing cold air 
Maintained room environment 
Thermocouple cold shoe 



TEMPERATURE 
100 F° 

70 F° 

75 F° 
291.5 K° 



TABLE 7 



For operation at low noise level, air velocity is selected 
to be 1800 feet per minute. 

The Reynold's number is determined from: 

N R£ = (4-6) 
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Using the air properties of Table 7 the corresponding 



value of H for each type of duct is calculated as: 



(n re>a = 


9285 


^RE^B 


8440 


^re'c 


8440 


* n re*d 


11973 



These values are indicative of turbulent flow. The calcu- 
lation of the convection heat transfer coefficients is 
based on the Dittus and Boelter equation: 



N nu = 0.0265 



(N )°‘ 8 (N ) 0,3 
RE pr 



(4-7) 



where and N^ r are the Nusselt's and Prandtl’s number 

correspondingly, given by: 




(4-8) 

(4-9) 



For a Prandtl's number almost equal to 0.7 the Nusselt's 
number, for the different type of air ducts, is: 
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35.1428 



(N NU* A 



°Wb 



(t Wc 



(n nu* d 



32.56 

32.56 

43.07 



The corresponding convection heat transfer coefficients are 
found from Equation (4-8) as: 



h A = 10.92 

h fi = 11.128 

h c = 11.128 

h D = 10.376 

Using the above values for the convection heat transfer 
coefficient and assuming aluminum as the fins material, with 
thermal conductivity 92 BTU/f t-hour-F° , the efficiency of 
the long and short fins, according to Equation (4-3), is: 

n,, = 0.9256 and n_ = 0.9716 

f 1 fs 

correspondingly . 

The rate of heat absorption is calculated from 
the following equation: 
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(4-10) 



V AT w^S s i 

where the subscript j corresponds to different types of air 
ducts of the heat exchanger. Analytically the exchange 
effective surfaces per thermocouple column are: 

S A = (5n fl + 3n fs + 1,08) 230 = 1983 - 2 cm2 = 2.135 ft 2 
S B = (5n fl + 0.26) 230 = 1124.24 cm 2 = 1.21 ft 2 
S C = 5n fi 230 = 1064.44 cm 2 = 1.146 ft 2 

S D = (5n fl + 3n fs + !- 232 ) 230 = 2018.2 cm 2 = 2.172 ft 2 

From the temperatures of Table 2 the logarithmic mean 
temperature of the cooling air, according to Equation (4-1) , 
is found to be: 



AT^= 15.0 



Substitution of the determined values for and hj in 
Equation (4-10) gives: 

q = (2.135-10.92 + 1.21-11.128 + 1.146-11.128 
+ 2.172-10.376) AT.. = 1081 BTU/hour 
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For the 12 thermocouple columns of the device, the capacity 
of the heat absorbing exchanger is found: 



q T0T = 12792 BTU/hour 

The above figure for the heat exchanger* s capacity offers 
a margin factor of 1.85 which ensures satisfactory operation 
under even, not considered, worse conditions. 

The air flow requirements can be determined from: 



C p (t in ' t out ) 



(4-11) 



for q = 12,792 BTU/hour, c p = 0.2 8 BTU/lbm-F° , 

t^ n = 100 F° , and t QUt = 70 F°, the rate of air flow is 

determined as: 



m = 1522.8 lbm/hour 

3 

For the air density f = 0.075 lbm/ft 

m = 5.64 f t 3 /sec 

3 . Hot Side Heat Exchanger 

This is also a finned type heat exchanger. It con- 
sists of a tube with radially introduced fins of two differ- 
ent lengths. The fin*s material is assumed to be aluminum 
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with thermal conductivity 92 BTU/hour-f t-F° . Following the 
same procedure as in the cold side heat exchanger the obtained 
results are displayed in Table 8. 



ITEM 


UNIT 


VALUE 


Required capacity 


BTU/hour 


13,503 


Hot shoe 


F° 


121.5 


Incoming air 

Temperature 


F° 


100 


Outgoing air 


F° 


110 


Log. mean difference 


F° 


16 


Tube diameter 


ft 


0.4462 


Long fin length 


ft 


0.1833 


Short fin length 


ft 


0.082 


Number of air ducts 


- 


24 


Duct equivalent diameter 


ft 


0.03784 


Duct length 


ft 


9.55 


Air velocity 


ft/sec 


86 


Reynold's number 


- 


20,860 


Nusselt's number 


- 


67.15 


Prandtl's number 


- 


0.7 


Convection coefficient 


BTU/ft 2 -hour-F° 


26.6 


Long fin efficiency 


- 


0.6296 


Short fin efficiency 


- 


0.882 


Calculated capacity 


BTU/hour 


18,983 


Security factor 


- 


1.4 


Volumetric rate of flow 


f t 3 /sec 


9.4 



TABLE 8 
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FIGURE 10 

ELEMENTARY T.C.-HOT SIDE 
HEAT EXCHANGER 
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V. THERMOELECTRIC GENERATOR 



A. ANALYSIS OF ELEMENTARY GENERATOR 

The elementary thermoelectric generator, which is operating 
under steady state conditions, consists of a thermocouple 
externally connected to a resistive load, as it is shown in 
Figure 11. 

The operation of the device can be described as follows: 
the hot junction of the thermocouple is assumed at temperature 
T c . The thermocouple arms are assumed to be adiabatically 
insulated, and heat exchange takes place only at the junctions. 
At the hot junction heat q Q is absorbed, and heat q^ is 
evolved at the cold junction per unit time. Under these 
conditions, a thermo-e .m. f . E = o( is established in the 
circuit and it maintains an electrical current I. The 
coefficient in the thermo-e .m. f . is the Seebeck coefficient 
of the couple 



where W is the available electrical power at the external 
load. 



*p(T) 4 c^Cr) | 



(5-1) 



From the law of energy conservation it follows : 




(5-2) 
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FIGURE II 

ELEMENTARY T.E. GENERATOR 
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w = I* Ro 
w - TE-I a R 
W = I*M-I e R. 



(5-3) 



where R is the total electrical resistance of the thermocouple. 

Equations (5-3) suggest the electrical equivalent of the 
thermocouple shown in Figure 12 . 



(K 



aaaaa 




Ro 



Figure 12: Electrical equivalent of thermoelectric 

couple . 



The general expression for the heat flow rate qg has the 
following form: 

v = TVCT.1I- MA f (r)£ „ , 1 -[ 
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This expression is reduced to a more practical form by 
introducing the estimate of the rate of heat flow in the 
conducting solid, as it is determined in Appendix A 



-}(T)Acr)iL 



dv 



- -\lfbll /\ 7. 



t 



.1 



Jlb + I 

(4t) r " 



s AT 



(5-5) 






•A 1 i SCO ( J ^CT)dT)dT 

6 = ~ 

AT f\ CUdT 

T c 

The plus or minus sign in equation ( 5— 5 ) allows for two 
possible directions of the electrical current: the minus 

sign representing the current in the P-type segment and the 
plus sign giving the current in the N-type segment. As it 
is shown in Appendix A, the value of £ is very close to half 
of the average resistivity for the range of temperature 

t h' t c* 

For convenience in further calculations, it is desirable 
to replace equation (5-5) with a less exact but much simpler 
expression : 




*c J 
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(5-7) 






T^T h 



where: 



rfj: ^ 



Combining equations (5-4) and (5-7) , the rate of heat 
flow at the hot shoe of the thermocouple becomes: 



v n(T>)I ' * ^ 

-I(VV> AT 



(5 



Introducing the Kelvin's relation, equation (5-8) becomes: 

V * Trt 1 + Bf * ] AT - k f [ fV* + ] - 1 A • - a *l at <5 

where 




and (5 



M)]olT 



- 8 ) 



“9) 



- 10 ) 
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The efficiency of the elementary thermoelectric generator 
is defined as: 



If I and AT are considered as constants, then the efficiency 
^ is a function of x and y which determine the geometry of 
the P and N type segments. 

Figure 13 shows the three dimensional representation of 
efficiency as it is interpreted in (5-11) , for a certain 
temperature range and as thermoelectric material. 

B. CONDITIONS FOR MAXIMUM EFFICIENCY 

Equation (5-11) is a function of geometric variables 
x, y and of the electrical current I, for a given temperature 
range. Using the Lagrange multiplier technique again, the 
necessary conditions for maximum efficiency are defined as: 




+ u = o 



Bx 5 r 



( 5-12a) 




(5-12b) 



^2L = o 



(5-12c) 
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where h is again the Lagrange multiplier and R = 
is the constraint function. 
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FIGURE 13 

THREE DIMENSION REPRESENTATION OF EFFICIENCY 
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Elimination of h from equations (5-12a) and (5-12b) gives 
the already found condition 



1 






\e,J 



Vz 



(5-13) 



Equation (5-11) can be written as: 



" 1 = 



ril AT -1 2 H 

^ T«1 + [-l p ,+ jkJ/n- |l 2 R.-lC« : r -4;)AT 



(5-14) 



Multiplying and dividing by R and introducing the condition 
of (5-13), the efficiency equation becomes: 



L *T„f - (Tat - i 

where 



(5-15) 



(5-16) 



Optimization of (5-15) with respect to gives the 
condition : 



* AT 



(5-17) 
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where : 



Q - ( i + -|r Tc( V 



-± _L \ V 2 

g 1 atJ 



Tcw= ±(T m + V) 



The optimum current is: 




u AT 

&•(<*+ q 



The optimum effective voltage is: 



(5-18) 



(5-19) 



V AT (5-20) 

When the conditions (5-13) and (5-17) are satisfied, the 
efficiency reaches its maximum value, which is found from 
equations (5-15) and (5-17): 






AT 

T* 



( G-il 






(5-21) 



The electrical power dissipated in the external circuit is 
found from equation (5-19) and (5-20) . 



y/- « a ftfrT) x 



(5-22) 



The formulae developed above show that the efficiency of 
a thermocouple and its power output increase with increase of 
the working range of temperatures. The thermocouple efficiency. 
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like the cooler C.O.P., depends only the the ratio T_/T„ 

C H 

and the quality factor G, which is slightly different from 
the quality factor D for the thermoelectric cooler. 

For the present application the working temperatures of 
the thermoelectric generator are varying as the temperature 
of the exhaust gasses varies for different power level 
operation of the automobile engine. In the present state of 
material development, no single thermoelectric material is 
satisfactory for use over the complete temperature range. 
Therefore, it is necessary to use different materials in each 
temperature range. One way to accomplish this, is to form 
each leg of the couple by joining different materials in 
series; that is, forming a sandwich structure of the appropriate 
thermoelectric materials for each temperature range. 

C. OPTIMIZATION OF THE SEGMENTED T.E . GENERATOR 

In Figure 14 two P-type materials, P^ and P 2 , are connected 
in series to form a positive thermoelectric arm and two N-type 
materials, N 1 and N 2 , form the negative arm of the thermocouple. 

The optimization problem is to maximize the thermal 
efficiency, subject to specification of the interface tempera- 
tures T and T over the equivalent length of thermocouple, 
x y 

L = In ( — ) , equivalent cross sectional area A p -(<4>l) , 

and the temperature averaged material properties. 

The reciprocal of the thermal efficiency, as defined in 
equation (5-11), is: 

4 q ThI + (Tm-T V) kpi* Oh"Ty) 1<vi ~ 5 1 0*pi 4 R/vi) - I ( g p“&A/)AT (5-23) 

"l ' l a R. 
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T c 




FIGURE 14 

ELEMENTARY SEGMENTED GENERATOR 
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Introducing the concept of equivalent length and cross section 
for the P-type and N-type segment as: 



u-uh 

-- U (f{) 

U, U%) 

L = + Q-?i - + in i 

A P = 

A/v= 



(5-24) 



The terms in the numerator of equation (5-23) have the 
following interpretations: 



Kp, = 



Kna - 
Rpi * 



— A N 

Ana n 

P ia 

^ A? 



Rri = 



P ^ 

U ' A* 



(5-25) 
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(5-25) 



£ = ± //< , i ( ' T) /A CT)cjTd 'r 

ri at /- Th 

1 l f1 C-r)dT 

*x 



^ At 



i / 6 "- C '° 1 CT)J TdT 

A -'TV T V 



r* 

L 



The electrical current I in equation (5-23) is given by: 



I 



E 

R + R 

o 



(5-26) 



where: R = + Rp 2 + ^1 + *N2 = *P + *n is ^ total 

electrical resistance of the thermocouple and E is the 
circuit voltage produced by the cumulative Seebeck effect 
in the positive and negative elements. The segmented elements 
may be considered as batteries connected in series. 

The voltage E is given by: 



E = <*' ( t h - V 



(5-27) 



where: 



cA' 



i 

AT 



[ / c V T)dT ^ c < p J t ) cJt ~J^ 



T* 

(T)cIt 
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Substituting equation (5-26) into (5-23), the latter becomes: 



i. _ ( T w -V) kn + (T h-Ty') 4 ^ Th RtEo ~ gcu-tRo^ ( s r>~ b-*Ko (5-28) 

V 



e‘R» 

OUR.? 



With the variables 6j> and 9,, defined by the equations: 




T w -T* 
T h - T c 



Dm - 



T» -Ty 
“T w -Tc 



(5-29) 



and multiplying the numerator and denominator of equation 
(5-28) by (R + R q ) 2 , equation (5-28) becomes: 

i . ie*(EM-nuO 

\ E a Ro (5-30) 

__ E ( g M 

E* Ro 

With the quantities R^, R^, K^, defined by the 

equations : 

V Mr 

t 

R^= R/vAn (5-31) 

4i = Jin 

Af 
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equation (5-30) becomes: 



/ / 

i__ (Tu-TjfGpU^Ap + ^W^Av] (~A^ h 



(5-32) 



+ ¥<,) (R?* 4 _ Cfp^ C 1 * ^ 

c < , C T W-T c . > ) Z ^o 



— - ( t H-~ t c) Lgp ^P + 6p ^?a Ryv ( Av) * 6*^*4 ( Ap ) * gyj^ 4 H if** gp) 

u” 



(5-33) 



+ * t h (*-*• ~^ 0 ) __ ( Rp4 +• £*0 _ ( g M - g >» a ) fi- 4 - ^ P ~) 

^'Cth-t c ) At«-V) 



Assuming / minimization of equation (5-33) 



B. " J x Ap 

with respect to 



, gives : 



From equation (5-34) the condition is obtained: 



A - 

a p ‘[e, \4< r’p \ 



Applying condition (5-35) to equation (5-33) , the latter 
becomes : 

i _ [ (©ptfpKpi) 

n " e 1 + 

(5-36) 

_ + ^ t h _ (Rph + Bam) _ AlAlAlL) 

°t' (.Th ~\) 2 R. 0 t< ' ^ t h - t c ) 

With the figure of merit defined as: 



AU 



U'fTa 



^[(OpR'p^Cg.CViV' 1 ] 

equation (5-36) becomes: 



(5-37) 



^ ’ 4 M Cl H --T t ) rt'(T H -T c ) 

(5-38) 

__ C Rpj 4 ^ (^P4 ~ 1 Cl* 1 

aRo 

1 c 

With the partial derivative of (— ) with respect to ~~ set 

L K. 

equal to zero, the following condition is obtained: 
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R. I ( TtL - Z-?2J±l ) \ 

L ^ A* X' ) 



Substitution of (5-39) into equation (5-38) gives the minimum 

_i_ or the maximum efficiency. 

\ 

Summarizing, the maximization problem of the segmented 
generator efficiency can be stated as follows: Given the 

values of ft ^ ^ ^y A T ij' T H ' T c' A p' and L determine 
the values of ^ j , R, R Q , I, and 0^ , where i = P, N 

and j = 1 , 2 . 

Solution of the problem requires that the following 
equations be satisfied: 

j _ d CTn-Te) 

Aw _ [ Bp kp< R.w 1 

M _ (*0 z t h 

4 



(5-40) 



(5-41) 
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Ho .! 


Tm( 


i 5 N 

j* - -spO I 


Vi 


(5-43) 


M 


(la 


(g P4 ” Gha "\ 






MM 








L = 


0P1 + 


p - p +P 

'-p2 ' '-M 1 t-Ux 




(5-44) 


R- 


R.pi + 


R-pa • + R/vi ■* 




(5-45) 



■^■1 R-P1+ 2.1 Rpa. ~ (^pa.- °S>i') "Th I " ( s pf* <Sp 2 )I +^ pi -|^(' t h-‘T^ = (5-46) 

= ^Tfz (T/ ‘ Tc) 

•|l R-A/x-t + f 6 AM + ■*■ ^' w 'l"^^ Ttt '' r '()(5-4 7) 

= Xr~T- ( T ^- T c) 

1 N-l' 

Equations (5-46) and (5-47) represent the energy balance for 
the interfaces between elements P^, P 2 and . They 

state that the amount of energy flowing to interface, plus 
or minus the amounts of energy generated or absorbed at the 
interface is equal to the amount of energy flowing from the 
interface. 

D. T.E. GENERATOR DESIGN 

The design of the segmented thermoelectric generator is 
implemented by computer aided iteration methods. Such methods 
have been coded for the IBM-360 as it is displayed in program 
Number one. The procedure consists of several steps: 
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1) 



2) 

3) 



For the selected materials, calculation of 
for i = P ,N and j = 1,2. 
Calculation of 0 p , 

Initial estimation of P f P 







from: 



(t ‘-^ 

^ il c t i - t =) 

L - + C** 



4) Calculation of R^ , R^, K pl , K p2 . 

r^I 

5) Calculation of — , from Equation (5-41) and for 

J? 

selected value of determination of 

6) Determination of figure of merit M from Equation 
(5-37) . 

R 

7) Calculation of — from Equation (5-42) . 

8) Calculation of R from equation: 

^ - ( fpi ^PA * + f 

9) Calculation of current I from Equation (5-40) . 

10) Calculation of both sides of Equations (5-46) and 

(5-47) . If their differences are greater than a 
predetermined figure, then a new set of 1 Pi' ■Hja 1S 
estimated and the procedure is repeated until the 
wanted accuracy is obtained. 
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11) Calculation of R^, K^. 

12) Calculation of efficiency from Equation (5-38) . 

The procedure is based on an "in advance" selection of 
the interface temperatures T x and with the criterion of 

safe behavior of the materials at the extreme operating 
temperatures . 

When the lengths of the P and N type segments are cal- 
culated for the extreme temperatures, the problem of re- 
estimation of the interface temperatures, for different hot 
junction temperatures, is encountered by modification of 
the iteration procedure in step 10. That is, for the calcu- 
lated lengths, and T^ are varying until Equations (5-46) 
and (5-47) are satisfied. 

The adverse operating conditions of the thermoelectric 
generator — high hot junction temperatures and extended 
operating range — impose restrictions in the selection of the 
best available thermoelectric material for energy conversion. 

All thermoelectric materials represent the best of 
their characteristics in a quite narrow range of temperatures. 
This problem is partially encountered by the segmented type 
generator design, where the one material acts as the comple- 
ment of the other in regions of poor characteristics. 

For the present application, the selected thermoelectric 
materials and some of their properties are displayed in Table 
9. 
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ITEM 


MATERIAL 


MAX 


OPE RAT 


. TEMPERATURE 


p i 


Ag Sb Te 2 






900 


(K°) 


P 2 


25% Bi 2 Te 3 


+ 75% 


Sb 2 Te 3 


600 


(K°) 


N 1 


Pb Te (25% 


Sn Te) 




900 


(K°) 


n 2 


75% Bi 2 Te 3 


+ 25% 


Bi 2 Se 


700 


(K° ) 



TABLE 9 



With proper dopants the maximum operating temperature of 
and materials can be considerably extended. 

The thermoelectric properties of the materials are 
displayed in Figures 15 through 21. The data have been 
taken from Reference [1] and [2] and the corresponding 
temperature relations developed using the "best fit poly- 
nomial" technique as it is shown in Program 2. 

With these data supplied, the computer runs through the 
iteration procedure described earlier, and solves for the 
optimum thermal efficiency and the optimum design parameters. 
The optimum solution is displayed in Table 10. 

The maximum and minimum values for efficiency, current, 
and voltage in Table 10 are associated with the highest and 
lowest hot junction temperature. 

Figures 22 through 24 shows the variation of efficiency, 
voltage, and current with respect to hot junction temperature 

T . 

H 



74 



D eg.~ ,L x 1 - Figure of Merit 
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FIGURE 15 
FIGURE OF MERIT 
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FIGURE 16 

THERMAL CONDUCTIVITY 
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FIGURE 17 

ELECTRICAL RESISTIVITY - OHM - CM x IO 3 
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Temperature - K° 

FIGURE 18 
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FIGURE 20 
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Temperature K° 



FIGURE 21 

ELECTRICAL RESISTIVITY - OHM - CM x 10 3 
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ITEM 



UNIT 



VALUE 



REMARKS 




-A P z (4>“O r 

V^P 

An* C^^/v 

R 0 / R 

R o 

R 

E 

I 

W 

n 



K° 

K° 



rad - cm. 

rad - cm. 

ohms 

ohms 

volts 

Amperes 

Watts 



min - 940 max 


Selected 


335 


Selected 


0.25845371 


Selected 


0.153 




0.10545 




0.22275 




0.0357 




o 

to 


Selected 


0.5 




0.1 




0.272846 




1.91 x 10“ 3 
7 x 10~ 3 





min 0.135 - max 0.225 
min 15.2 - max 23.6 

min 0.44 - max 1.41 

min 6% - max 9% 



TABLE 10 
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FIGURE 22 

CURRENT (AMPERES) VERSUS HOT JUNCTION TEMPERATURE 
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FIGURE 23 

EFFICIENCY VERSUS HOT JUNCTION TEMPERATURE 
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Temperature K° 

FIGURE 24 

VOLTAGE VERSUS HOT JUNCTION TEMPERATURE 
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E. ESTIMATION OF AVAILABLE HEAT POWER FROM THE EXHAUST GASES 

The recently introduced catalytic converter in the 
exhaust system of automobiles, offers the advantages of 
additional release of heat and higher temperatures, by 
oxidation of the unburned CO and H 2 - 

Heat, also, can be extracted by cooling down the gasses 
which come out of the output flange of the converter. 

The thermodynamic analysis which follows is based on 

3 

data provided by FORD motor company. A 2300 cm engine is 
examined and the results of the exhaust gas analysis before 
and after catalytic converter are summarized, for three 
different power levels, in Table 11. 

The analysis is divided into three parts with the neces- 
sary logical sequence. 

The first part contains preliminary calculations, as 
the moles/min and pressure of exhaust gasses. The second is 
involved with the calculation of the available heat energy 
from the catalytic action of the converter, and the third 
is referred to the estimation of the possible extracted 
heat energy by cooling down the gasses after the catalytic 
converter. 

1. Preliminary Calculations 

a) From Table 11 at idle speed, the engine produces 

3 3 

exhaust gasses at a rate of 10 ft /min or 0.2832 m /min. 

One mole of any gas at the state of 1 atm and 273 K° obtains 
a volume of 0.0224 m 3 and at standard state — 1 atm, 298 K° — 
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0.02445 m^. Therefore, the flow of the exhaust gasses in 
terms of gram-moles per minute is 11.583. 

Based on the percent mole analysis of Table 11, 
the actual amount of the several constituents of the exhaust 
gas can be summarized as follows: 





CURB IDLE 


60 MPH ROAD LOAD 


PEAK POWER 


GAS 


gram-mole/min . 


gram-mole/min . 


gram-mole/min 


co 2 


0.80269 


8.2700 


17.316 


CO 


0.4981 


0.3475 


12.787 


h 2 0 


1.0500 


8.548 


24.509 


H 2 


0.1865 


o 

o 


5.238 


°2 


0.64864 


3.3358 


7.992 


N 2 


8.39707 


48.995 


198.472 


TOTAL 


11.583 


69.497 


266.404 



TABLE 12 



b) The partial pressures of the gasses can be 
calculated from: 

r " V Tot (5-48) 



where : 



= gram-mole of the gas 
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R = gas constant 
T = gas temperature 

V TOT = vol 111116 occupied by all gasses, comprising 

the mixture at temperature T. 

With taken from Table 12 the partial pressures 
were calculated for three different power levels, and the 
results are shown in tabular form as follows: 





CURB IDLE 


60 MPH ROAD LOAD 


PEAK POWER 


GAS 


atm 


atm 


atm 


co 2 


0.06931 


0.1190122 


0.066007 


CO 


0.04301 


0.005001 


0.048 


h 2 o 


0.09067 


0.123013 


0.09201 


H 2 


0.0161 


0.0 


0.1967 


°2 


0.056007 


0.048005 


0.030003 


N 2 


0.725045 


0.70508 


0.7451 


TOTAL 


1.000142 


1.0001112 


1.00079 



TABLE 13 

Based on the data of Table 13, it is reasonable to 
consider that the chemical process in the catalytic converter 
takes place under constant pressure conditions. 
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2 . Released Heat From Converter Action 



The change in enthalpy produced by the procedure 
of a chemical reaction from reactants in their natural state 
at 25° C and 1 atm to products in their natural state at 
25° C and 1 atm is referred to as the "standard heat of 
reaction" at 25° C. 

The "standard heat of formation" is the change in 
enthalpy resulting from the formation of a compound at 
25° C and 1 atm from the elements in their natural states 
at 25° C and 1 atm. The "standard heat of combustion" is 
the change in enthalpy resulting from the reaction of a 
compound with elemental oxygen each initially at 25° C and 
1 atm to produce specifically defined products at 25° C 
and 1 atm. 

The heat of a given reaction may be obtained by 
constructing a set of formation and combustion reactions 
such that the sum of the individual reactions is the reaction 
desired. The sum of the enthalpy changes of the individual 
reactions will then be the standard enthalpy change of the 
given reaction. If the heats of formation data are available 
for all compounds involved, the standard heat of reaction 
may be obtained by subtracting the sum of the heats of for- 
mation of the reactants from the sum of the heats of forma- 
tion of the products. Normally, heats of formation are 
available for inorganic compounds and heats of combustion 
are available for organic compounds. 
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Since the initial and final states differ from 25° C 



and 1 atm, the change in enthalpy may be evaluated by a 
series of steps so chosen that the enthalpy change of each 
individual step can determined. These steps are shown on 
the following diagram. 

The change in enthalpy for steps 1 through 3 and 
5 through 16 are calculated from: 



-- Hij r\. 



dT 



(5-49) 



where M^ is gram-mole and c^ is the molal heat capacity of 
the gas. The molal heat capacity is a strong function of 
temperature and can be expressed as : 

C PL = 1 +c i. T (5-50) 



where T is in K° . 

The coefficients o^, b i and are displayed in 
Table 14. 



COMPOUND 


FORMULA 


0(i 


b ± x 10 3 


c^ x 10 


Carbon dioxide 


CO 2 


6.214 


10.396 


- 3.545 


Carbon monoxide 


CO 


6.420 


1.665 


- 0.196 


Water 


H ? 0 


7.256 


2.298 


0.283 


Oxygen 


o 2 


6.0954 


3.2533 


- 1.0171 


Hydrogen 


H 2 


6.9469 


- 0.1999 


0.4808 


Nitrogen 


"2 


6.4492 


1.4125 


- 0.0807 



TABLE 14 
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Step 4 is calculated from tables where the heat of 
formation for different compounds is given. The change in 
enthalpy is obtained by subtracting the heat of formation 
of the reactants from that of the products . 

Table 15 displays the change of enthalpy for each 
step, for three power conditions: 





CURB IDLE 


60 MPH ROAD LOAD 


PEAK POWER 


ENTHALPY 

CHANGE 


Calor/min 


Calor/min 


Calor/min 


H 1 


- 1340.11 


- 1418.0 


- 48025.0 


H 2 


- 955.36 


- 738.7 


- 36487.0 


H 3 


- 492.60 


o 

o 


- 20702.0 


H 4 


-46303 .00 


-21285.0 


-824828.0 


H 5 


+ 2394.40 


+ 2114.0 


+ 70449.0 


H 6 


+ 728.60 


o 

o 


+ 25102.0 


H 7-12 


+ 741.00 


- 930.0 


- 6788.0 


H 8-13 


+ 741.30 


- 1157.0 


- 6354.0 


H 9-14 


+ 4859.00 


- 4129.0 


- 44372.0 


H 10-15 


+ 714.00 


- 282.0 


0.0 


H ll-16 


+ 3.40 


0.0 


+ 885.0 


TOTAL 


-38909.37 


-27826.0 


-892886.0 



TABLE 15 
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3. Estimation of Available Heat from Hot Gasses 



The gasses flowing out of the catalytic converter 
have a high thermal content, which partially can be extracted, 
introducing a heat exhanger in series with the catalytic 
converter. 

A possible arrangement of converter with heat 
exhanger in series is shown in Figure 25. The change in 
enthalpy, because of the cooling process of the exhaust 
gasses, can be estimated from equation (5-49) . 

For a temperature drop 200 F° , the change in enthalpies, 
for different load conditions are shown in Table 16. 





CURB IDLE 


60 MPH ROAD LOAD 


PEAK POWER 


ENTHALPY 








DROP 


BTU/min 


BTU/min 


BTU/min 


H co 2 


1737 


10119 


31173 


H co 


5 


15 


3468 


H °2 


102 


2400 


0 


h h 2 o 


1355 


7846 


30700 


H « 2 


0 


0 


1063 


% 


6894 


35080 


169892 


TOTAL 


10090 


55460 


236296 



TABLE 16 
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FIGURE 25 

CATALYTIC CONVERTER AND HEAT EXCHANGER COMBINATION 
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In summary, the total amount of available heat power 
from the exhaust gasses of the examined engine, under the 
assumed conditions is determined as: 

1) Idle condition 38909 + 10090 = 50000 cal/min 

= 3486 watts 

2) 60 MPH ROAD LOAD 278264 + 55460 = 83286 cal/min 

= 5808 watts 

3) PEAK POWER 892886 + 236296 = 1129182 cal/min 

= 78741 watts 

F. FINAL CONSIDERATIONS ON GENERATOR 

The number of required thermoelectric couples for a 
typical design of the generator is based on the available 
thermal power which corresponds to 60 mph road load of the 
automobile engine. 

The design specifications of the thermoelectric generator 
are shown in Table 17. 
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ITEM 


UNIT 


VALUE 


REMARKS 


Available thermal power 


watts 


5808 




t h 


K° 


900 


SELECTED 


T c 


K° 


335 


SELECTED 


Efficiency 


- 


8.9% 




Electrical output power 


watts 


517 




Current 


amps 


22.3 




Voltage 


volts 


23.2 




Number of couples 


- 


43.2 




R i 


cm. 


8.00 




R 2 


cm. 


10.36 




4>p 


rad 


0.2 




cf’N 


rad 


0.1 




-z 


cm. 


1.00 




cj?insulation 


rad 


0.1 





TABLE 17 

\ 

The maximum working temperature of the selected thermo- 
electric materials is restricted to 900 K°. It is possible 
to have higher temperatures in the wall of the catalytic 
converter when the engine is overloaded. In order to protect 
the hot junction thermoelectric materials from damage, a 
protective system must be provided. Such system could be a 
single by-pass line for the exhaust gasses, automatically 
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operating when a certain temperature is achieved, or intro- 
duction of a heat resisting material between the converter's 
wall and the hot junctions of thermocouples, which will 
insure a sufficient temperature drop. 

The cold junction temperature was arbitrarily selected 
to be 335 K°. It can be even lower, depending on the heat 
rejection system which is employed, i.e. water cooling, forced 
or free air convection. 

The hot junction temperature is varying according to the 
developed power level of the engine. This results in the 
variation of all design specifications of the thermoelectric 
generator, as delivered useful power, voltage, current, etc. 
For this reason, an electrical power conditioning device 
becomes necessary. The function of this device depends on 
the special application of the generator. That is, if the 
generator is employed as the electrical power source of the 
automobile, constant voltage and variable current is required, 
and when the generator supplies part of the electrical power 
to the thermoelectric cooler, constant current and variable 
voltage is desired. 

Contact resistance and conductor losses were not taken 
into account, because the estimation of the former is a matter 
of manufacturing procedure and the latter depends on the 
proper material for the applied contact technique. 

Nevertheless, a reasonable figure of contact resistance 

2 

is of the order of 10 micro-ohms per cm , which is very small 
compared to the resistivity of the thermoelectric materials. 
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Figure 26 shows the arrangement of catalytic converter with 
the attached thermoelectric generator. 
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FIGURE 26 

T.E. GENERATOR ASSEMBLY 



100 



VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 



The obtained figures for the thermoelectric cooler and 
generator could be considered as quite promising at least 
from the standpoint of feasibility of the application. 

The thermoelectric cooler for 7000 BTU/hour was found 
to require 1953.0 watts of electrical power. From this figure 
50 watts are dissipated for the required air blowers. If the 
four units, which comprise the device, are fed in parallel, 
then an external conventional generator, coupled to the 
engine, is needed with 33. volts and 57.6 amperes as rating 
characteristics . 

When the thermoelectric generator and the cooler are 
considered as one system, then the required external electrical 
power is reduced by 18 to 25%, depending on the operating 
power level of the engine, and correspondingly the coefficient 
of performance of the cooler is increasing. 

A direct comparison between the conventional air condi- 
tioning system and the thermoelectric cooler-generator scheme 
is not easy, because there are not available enough data for 
the latter. In Table 18 an attempt for comparison is made in 
the most broad characteristics. 
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THERMOELECTRIC 

SYSTEM 



CHARACTERISTIC 



CONVENTIONAL 

SYSTEM 



Volume 

Weight 

Complexity 

Maintenance 

Noise level 

COP 

Working fluids 



Both systems have almost the same 
worst best 

complex simple 

required almost none 

relatively noisy quiet 

Both systems are in the same level 
yes no 



TABLE 18 



From the standpoint of manufacturing cost, the existing 
figures show that the thermoelectric system is two to three 
times more expensive from the corresponding conventional. 

It must be mentioned that the cost of the conventional system 
is referred to the "series production cost" and the corre- 
sponding for the thermoelectric system, to experimental 
constructions. It is expected that whenever the thermoelectric 
system will be manufactured in series, the cost will decrease 
enough and eventually it will be compatible or even cheaper 
from the conventional system, since the thermoelectric material 
itself is not expensive. 

The analysis showed that the efficiency of the generator 
depends strongly on the thermoelectric material. The obtained 
figures — 6 % minimum and 9% maximum — for the generator's 
efficiency, are considered satisfactory compared to 
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efficiencies reported already in the literature. Bates and 
Weinstein developed a similar design using PbTe and SiGe 
which complemented each other over the temperature interval 
from 50C° to 1000 C° . Their generator achieved an efficiency 
of 8.7 percent with a hot junction temperature of 800 C° 

(1073 K°) . The present design shows an efficiency of 9% 
with T h = 940 K° and T c = 335 K°. 

The requirements for the primary electrical source of a 
small automobile vehicle are satisfactorily fulfilled from 
the specifications of the segmented type thermoelectric 
generator. 

Table 19 displays the electrical power requirements for 
a small size car. 



ITEM 




WATTAGE 

(12-Volts) 




high beam 


45 


Head lamps 


low beam 


40 


Front lamps 




21 


Tail lamps 




21 


Other lamps 




10 


Wiper motor 




10 


Electrofan motor 


11.5 


TOTAL 




158.5 



TABLE 19 



The thermoelectric generator provides 230 watts in the 
idle condition of the engine and 517 watts in the 60 mph road 
load. So, from this simple examination, the thermoelectric 
generator can replace the conventional alternator, offering 
the great advantages of fuel economy and elimination of 
moving parts. 

The realization of the discussed applications requires 
further theoretical and practical work to be done. The 
systems must be examined thoroughly from the following 
standpoints : 

1) Stress analysis of T.E. couples for T.E. generator 
and T.E. cooler. 

2) Thermal and electrical insulation. 

3) Contact techniques . 

4) Fabrication of thermoelectric pellets. 

5) Power conditioning circuits for T.E. generator. 

6) Control system for the automatic co-operation of cooler, 
T.E. generator and conventional alternator. 

It must be stressed that the present work was developed 
under the scheme of the feasibility study instead of a 
detailed design. 
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APPENDIX A 



APPROXIMATE ANALYTIC SOLUTION OF HEAT BALANCE EQUATION 



The heat balance equation is: 

VT)m— iIs(T)il +[ 2 JH) = 0 (a-u 

4 dc x dT \dr-j ^ J dr dr ^nr 

where: }(,Y) #■ f(T) r €(t) r are functions of temperature. 

Equation (A-l) can be written also as: 



"^(J) 'tor'^A-I -t~ rii 1 IsC^-wril 4 I 2 p( T ).-0( A -2) 

dr" 1 dT x dr 1 dr dr ^ 



dT dT dT 

Expression V -r- can be written as — 7 - or — 

dr d(W) 



If the variable x is introduced such that x = In r then: 



r 



_ c!t 

dr d * 



(A-3) 



dT ± dr 
dr r dy 



(A-4) 



J rdT 1 _ d_r_L dT 1- +_L d^T d*. 

dr Ldr ] " dv L r d> J dr *-r] dx f dn l dr 



(A-5) 
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(A-6) 



i_ dr + _J_ <J*T 

d ' T'dx r 1 d x z 



.-•* <f T „ J*t _ 

dr 1 dx 1 ~ dx 



(A-7) 



Applying transformation (A-7) to equation (A-2) , it becomes: 

+ f e«=° 



or: 



)l T W ~ + IsCTiOfliL +1 7 pCt) = 0 

dX T dl Mx/ dx C 



Dividing by the product ml equation (A- 8) becomes: 

A j_ + JCt\Ct)3 2 ^ (dj ^ + eQ0! h ^ j p(-r 3 ^ 

X dx z dT 2 M*/ d* on_ “ 



( A-9 ) 



The following transformation is introduced: 



-H = - 



^CT^nnrL 

X 



dT 

dX 



(A- 10 ) 
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Differentiation with respect to T, equation (A-10) becomes: 



_ _ ttl_ d^CrV) 
cil* “ 1 dT 




I ^ dT 



or 




dT 'd* ' I dT ' d-x j j 




(A-ll) 



Applying the transformation (A-ll) to equation (A-9) , the 
latter becomes: 



Through the attempted transformation the order of the original 
differential equation was reduced but the non-linearity 
remained, so the solution of it at least is not obvious. From 
this point, the analysis of the problem is based on the 
assumption that in a first approximation the electric effects 
are independent. That means, it is assumed that the Thomson 
effect is not interacting with the Joule effect. 

From this standpoint, the electric effects can be studied 
separately and their influence can be added finally in the 
interpretation of the formulae for the rate of heat flow in 
the two ends of the semiconductor element. 



dT 



ch + ,ft T)XT) . sCT -) = Q 
dT 2. 



(A-12) 
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1. Influence of Thomson Effect on Heat Conduction 



It is assumed that the Joule effect is not present in 
the system. Then it can be set p(T) =0. 

Equation (A-12) becomes, for one possible current direction: 

ifL. 4 6cn=o (a-13) 

dT 

Integrating and introducing the transformation (A-10) , 
equation (A-13) becomes: 

= "B- l/ S'CT'jdT (a-14) 

**c. 

where B is the integration constant which must be determined. 
Transforming the independent variable x to the original r, 
equation (A-14) becomes: 

- - b- ij so) <Jt (a-is) 

The boundary conditions of the system are: 




Under these conditions, the rate of heat flow in the two ends 
of the element is given by: 
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(A-16) 



-V r ^ W ' r T I | „ « B-l/ 6 -CT) Jt 

6 t ir= 



Equation (A-15) can also be written as: 



~')lCO‘ irn ' r — = 'B (l- — [ 

V> Au * 



(A-17) 



Integrating both sides: 



B = ' W1 ; 


r ^C-r)dT 






L 1 -%J 


f T 6CT)dT^j 



(A-18) 



The following notations are used: 



Q 

q T 



rate of conduction heat flow = 




rate of Thomson heat flow = 





From equations (A-16) it can be seen that B^Q . 

The equality sign holds for the low temperature end. 

Also, for normal conditions, Q> | q | and consequently, 

l q Tl y 

— - — < 1 , which justifies the required conditions in order 

D 
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to expand the fraction 



1 



in a series of the 






following form: 






Equation (A-18) now can be written as: 






Integrating term by term and after some algebra, equation 
(A-20) can be written as: 



B ^ Q [ 1 + 





(A-21) 



where : 




t) ( 4 0 6 Ct)cIt) ol-r 

Ik T)dT 

% 



Vc 



for i = 1,2, ... ,n 



From equation (A-20) it can be set, for a first approximation, 
B = Q. Then, B for the second approximation becomes: 
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It is sufficient to keep only the first two terms of equation 
(A-22) , since higher order terms are quite small. 



B = fi> * Q 



1~T 

/ T « 

J ^CT)dj 

X. 



(A-23) 



Substitution of (A-23) in equation (A-16) yields: 
















yt. /' Ta i r J ecr > L 

— / - J -±5. 

/tUN L 'I X /T* 



«■(%) 



r H 

/ T vt)j t 



OD 



»«.' M10 



T “ / 3CO /cCOd-Tch- 

-v J r 4 - 



(A-24) 






2 . Influence of the Joule Effect on Heat Conduction 
Here it is assumed that the Thomson effect is not 
present, so setting (5 (T) =0 equation (A-12) becomes: 



CIt)^Ct) 
dT ^ 



= o 



(A-25) 



Rearranging and integrating, it is obtained: 



D 1 — ay^CT^ COdT 

where D is the integration constant. 

Introduction of the transformations: 

dT 

2 

I d* 




in the above equation yields: 

fV 1 ^ C ' T)dT ] 1fa (A " 26) 

The appearance of two signs in equation (A-26) leads to the 
conclusion that the temperature distribution along the element 
will present a local maximum T m at a point corresponding to 
some radial distance R^. 

Rearranging equation (A-26) , integrating both sides 
and applying the boundary conditions, the following expression 
for D^is obtained: 
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and for the point (R , T ) 
m m 




(A-2 8) 



There is a certain value of the current which could 



establish the condition of no back heat flow in the hot end 
of the element. That means, it is possible to have the point 
(R m , T m ) shifted to coincide with the point (F^ T fl ) , applying 
a certain value of current. 

This value of current can be found from equations 
(A-2 7) and (A-2 8) and it is: 



The above result might be important for the thermo- 
electric cooler, which should operate under the most efficient 
conditions. 

The rate of heat flow at the two element ends is given 
by the following equations. 





( A- 3 0 ) 
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It is obvious from equation (A-26) , in order for the 
left side to be real, the following condition should exist: 

D* > 2.X 1 f ftT") }CT) dT (A-31) 

for the range of temperature under consideration. 

Also, from equation (A-26) by integration of both sides 



on f 


, T M _ 

/ XCT)dT j 


(A-32 ) 


M-» 1 


*c- 





The condition declared by equation (A-31) justifies 
the series expansion of the integrand of equation (A-32) , in 
the binomial form. 



( i - j pCTX^CO dnrj ^ - C' l4y ) 

where p =-j- and x = - ~j^r j' ^" T ^C-'T) dT 

1 

( 1 + x) P = i Px + KP-Qx 1 ^ P(p- 0 (p- 2 ~)x s . + 

2. 1 "i> l 



(A-33) 



'n.l 

i . Xv 4 _S_ JL* , _15_ *?-*...+ X J<? 

1 Z H 2-! * 5i a " 1 nil 



114 



Application of the above expansion in equation (A-32) and 
integration of it term by term would yield, after algebraic 
manipulation, the following expression for D: 



Q n It 2i c* 



5L0. 







(A- 34) 



where : 



n I ' 



q 



Jn 







±_ 

TL 



Following the same procedure as before, one can set 
D^= Q as a first approximation. Applying this value of D L in 
equation (A-34) , can be calculated in a second approximation 
as : 



D = Q [ l + -1^- ^ Jaa. _L 4 . . 1 



(A-35) 



The higher order correction in equation (A-35) might be 
considered to be quite small to be taken into account in an 
approximation solution. 
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(A-36) 



It is sufficient to approximate D ( by: 

Applying the above expression for in equations (A-30) the 
rate of heat flow at the two ends of the element, now it is 
given by: 






(A-37) 






<h 



- (/v^r 



3. Composite Solution 

Based on the assumption that the electric effects — 
Joule and Thomson — are non-interacting, superposition of 
the obtained results yield the following composite solution 
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for the rate of flow of conduction heat in the semiconducting 
pellet: 

™ U J A Tm t CJ^Crtch) 

, / 6-Ct) f ~\Cj~)ATdT 
. T Ac A. (A _ 38) 



-}C0 



mn'T 



Jr 

dr 






/it rt 

A ^eodT 

,TH ^-T 

yn r i t z Q (—*) AW 
= r 2 - J 2 — ^yh fi -AUJ 



(X3^ t ) 

T 

/TM 

X_ ^COdT 



(A-39) 



Equation (A-38) is quite complicated to be used in the analysis 
of the thermoelectric cooler or generator, and it can be 
replaced by a simpler but less exact expression as: 



_Y T) ^ r ll | = fv-odT _ ± I 7 p ilS _ l£ AT 

^ dr [r=^ ^ ' w 



(A-40) 



where: 



__l X 6CT) /3 CT)dTjT 

^ ^CT)dT 
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The introduced error in equation (A-40) has been 
investigated using computer techniques and was found to be in 
the order of +5.8% and -3% for the resistivity terms of (A-40) 
for P and N type semiconducting pellets. 

Analytically, a semiconducting SiGe allow has been 
used with electrical resistivity and thermal conductivity as 
shown in Figures 27 and 28. For T R = 645 K° and T c = 400 K° 
the following results were obtained: 



These small errors have a very slight effect in the 
calculated amount of the conduction heat rate of flow as it 
is calculated by (A-38) and (A-40) . The calculated error 
was of the order of 0.2%. 

The validity of the assumption of non-interacting 
electrical phenomena, in a first approximation, can be 
justified by replacing the arbitrary relation C“0 

with constant values in equations (A-38) and (A-39) . Such a 



P-TYPE 



N-TYPE 




£ = 1.12535 x 10 -3 £ = 0.64771 x 10 

f = 1.275 x 10" 3 j f = 0.61 x 10 -3 




where : 
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Temperature K° 

FIGURE 27 

THERJ.1AL CONDUCTIVITY x IO 3 - WATT/CM/DEG. 
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Temperature K° 

FIGURE 28 

ELECTRICAL RESISTIVITY - OHM - CM x IO 3 



120 



substitution gives the exact form of the rigorous solution 
of the simplified differential equation as obtained in 
Appendix B. 
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APPENDIX B 



SOLUTION OF SIMPLIFIED HEAT BALANCE EQUATION 



The heat balance equation in the thermoelectric element 
with constant properties was defined in Section II as: 



ill + _L il 
dr 1 r dr 



where B 



fi* 

Ti(4xf 




Multiplying by r 



2 



equation (B-l) becomes: 



(B-l) 



r-\£l + r-iL 4. B = o (b- 2) 

dr' 1 - dr- 

Introducing the transformation x = In r , equation (B-2) 
becomes : 



CpI - _ B (B-3) 

d>* 



Integration of equation (B-3) yields: 



_ = -Bx -V (B-4) 

d* 



T= C 2 4C,X+iBP — - (B-5) 
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Replacing x with In r equation (B-5) becomes: 



T- C 2 -v C'vlV* 1 & ( 


(B-6) 


The initial conditions for the system are: 




r = R 1 T = T r 




O 

II 

E-* 

<N 

II 

u 




Introducing the above conditions to equation 


(B-6) , it yields: 


T w 


(B-7) 


Tc-c^c.k 


(B-8) 


From equations (B-7) and (B-8) the constants 


of integration 


and C 2 are determined as : 




^ At 


— (B-9) 




C z =Th-C a + 1| [CCO]* 


(B-10 ) 
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The rates of heat flow at the hot junction (T = T„) and at 

n 

the cold junction are defined as: 



\ * lf 4x 



v Lc^ 4t - tI,f ^ 



(B-ll) 



(B-12) 



The temperature distribution along a SiGe segment is shown 
in Figure 29. 
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length - cm 



FIGURE 29 

TEMPERATURE DISTRIBUTION IN A SiGe SEGMENT 



PROGRAM NO. 1 

THIS PROGRAM CALCULATES THE CHARACTERISTICS OF THE 
SEGMENTED THERMOELECTRIC GENERATOR 
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CALCULATION OF CONDUCTIVITY INTEGRALS 

LP1 = 60P*DTX1UB1 P*(DTX12/2. I *8 2P *DTX 1 3/ 3 
LP2=BPO*DTX21 *BP 1 * ( 0TX22/2 . ) +B P2 *DTX23/ 3 
LN1 = bON*DTY1UBIN*(DTY12/2. )*-B2N*DTY13/3 
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>BC = G P+ GN 

AV( I )= (TH( I i+TC( I ) i/2 
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